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>v Abstract 
o3 



The Hamiltonian dynamics of spherically symmetric massive thin shells in 
the general relativity is studied. Two different constraint dynamical systems 
representing this dynamics have been described recently; the relation of these 



, two systems is investigated. The symmetry groups of both systems are found. 

^ . 

New variables are used, which among other things simplify the complicated 

00 
0\ 



system a great deal. The systems are reduced to presymplectic manifolds T\ 
and T2, lest non-physical aspects like gauge fixings or embeddings in extended 
phase spaces complicate the line of reasoning. The following facts are shown. 



Ti is three- and T2 is five-dimensional; the description of the shell dynamics 

H : 

by Ti is incomplete so that some measurable properties of the shell cannot be 

predicted. T\ is locally equivalent to a subsystem of T2 and the corresponding 

local morphisms are not unique, due to the large symmetry group of r^. Some 

consequences for the recent extensions of the quantum shell dynamics through 

the singularity are discussed. 
04.60.Ds, 04.20.Fy 
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I. INTRODUCTION 



The dominating physical problem of the theory of gravitation is the gravitational collapse 
and the inevitable singularity. That is the point, where the classical theory breaks down; 
one expects that the quantum theory will help. 

The project of which the present paper is a part focuses on simplified models of grav- 
itating systems, whose quantum mechnics can be constructed without much technical and 
conceptual difficulty We hope that such models can help us to find quantization methods 
that are, in short, a) gauge (reparametrization) invariant and b) liberated from semiclassi- 
cal thinking. Indeed, the gauge invariance is an issue, because a gauge fixation leads much 
more easily to gauge- dependent results in quantum gravity than, say, in the Yang-Mills field 
theory. This has to do with the so-called problem of time 0. The term 'semi-classical 
thinking' is to express uneasiness about the extensions to the gravidynamics of the naive 
assumption that every dynamics is a theory of motion of some objects in some spacetime (cf. 
'covariant quantization' ||, string theory [|J, or 'effective field theory' ||). As the space- 
time is itself dynamical in the presence of gravitation, the only method complying with the 
assumption seems to be an expansion around a classical solution, and so it must be, it seems 
at least, a sort of WKB approximation. 

Our pet model is the spherically symmetric thin shell of dust in the general relativity. In 
refs. f|, [0, || and |§, a sufficiently simple super-Hamiltonian for this system was guessed 
from the equations of motion, so that the model could be readily quantized. We shall call 
this method the Warsaw approach^. This quantum theory lead to some unexpected results: 

1. Existence of a unitary scattering theory; after a contracting phase, the shell goes 
through a more or less probable intermediate state of black and white hole and then 
spacetime wave packet expands again; the probability of the hole stage depends on the 



1 A systematic and general exposition of this method (Ref. ||) was first given at a Banach Center 
Workshop, Warsaw. 
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energy of the shell. 

2. This unitary regime exists for all values of energy but only if the total rest mass does 
not exceed the Planck mass; what happens with more massive shells remains unclear. 

3. There is no spacetime geometry that can be associated with any particular scattering 
process: the wave packets contain linear combinations of wave functions that describe 
spacetimes with contracting shells and black holes as well as spacetimes with expanding 
shells and a white holes. 

The method, however, has several weak points. 

1. The way from a dynamical equation to a Hamiltonian principle is non-unique and it is 
unclear how much the final quantum results depend on a particular guess. Some kind 
of uniqueness has been shown in Ref. ||, but this is not yet completely reassuring, as 
we shall also see in the present paper. 

2. The great simplicity of the super-Hamiltonian is achieved by a very special choice of 
variables: they are the coordinates of the shell in the subspacetime at one side of 
the shell. This subspacetime of description (left subspacetime) lies on the opposite 
side of the shell than the subspacetime where the observers spend most of their time 
(right subspacetime). The drawback is that some observable properties of the shell 
are not contained among (or calculable from) these variables. A controversy arises as 
to whether and where the observers will see the expanding shell. Undetermined are 
for example the scattering time delays of the expanding packet with respect to the 
contracting one (usually given by the derivatives with respect to energy of the phase 
shifts of the S-matrix), which should in principle be measurable. The possibility was 
even discussed that the wave packet expanded into a different right subspacetime than 
where it originally contracted so that the scattering time delays made no sense. 

3. The method of self-adjoint extension of Hamiltonian operator was applied in Refs. ||, 
[0], and || to make the quantum evolution complete (unitary). This method seemed 
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to some physicists too formal and suspicious to cope satisfactorily with the problem 
of the singularity. 

The present paper is an attempt to deal with all three problems. The analysis will be 
based on Refs. |TIJ and ]TI]; we shall refer to this approach as the Potsdam oneQ. The 



Potsdam Hamiltonian action for spherically symmetric massive shells has been derived from 
first principles, so one can say something about the first problem. The Postdam dynamics of 
the shell is described by its coordinates in both subspacetimes left and right from the shell, 
and so it contains complete information about the motion. This enables us to say something 
about the second problem. Finally, we shall find that the method of self-adjoint extension 
is incomplete in certain respect, and so there is some information about the third problem, 
too. There is, of course, a strong motivation for the Warsaw approach: this is the etremely 
simple form of its action. 

The plan of the paper is as follows. In Sec. [H], we shall briefly collect the results of 
Refs. and [[LI]] that will be needed, so that the paper becomes self-contained. New is 



the systematic use of the Kruskal coordinates and more or less complete description of the 
symmetries that the systems admit. Sec. flTJ is based on the idea that the 'naked' or minimal 
mathematical structure that underlies a given constrained system is the constraint manifold 
together with the induced presymplectic form. According to this idea, two constrained 
systems are equivalent, if they define the same presymplectic manifolds. This notion is 
completely gauge (and reparametrization) invariant; it is moreover independent of the details 
of imbedding of the constraint manifold in an extended phase space, which is not physical and 
not unique. In Sec. [Ill], we reduce the Warsaw and the Potsdam descriptions to these 'naked' 
forms. The Warsaw presymplectic manifold (TijGi) is three-dimensional and the Potsdam 
presymplectic manifold (T2, O2) is five-dimensional. As a byproduct of the transformations 
made in this section, 'miraculous' new variables are found which considerably simplify the 



2 Ref. lol was written down at the Albert Einstein Institute, Potsdam. 
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Potsdam action. 

In Sec. [IV], we pick up a three-dimensional subsystem (Te, ©#) from the Potsdam system 
(T2, 62) that has a chance to be equivalent to a particular Warsaw system. Concretely, the 
Schwarzschild mass of the left subspacetime must be fixed, and the corresponding cyclic 
coordinate must be excluded. Sec. [VJis devoted to the search of a morphism of presymplectic 
manifolds that would map (T\, ©1) onto (r^, ®e)- We set up a partial differential equation 
for this map and solve it in suitable coordinates. A careful study of the results reveals that 

1. the two presymplectic manifolds are locally but not globally equivalent, and 

2. the local equivalence maps (morphisms of presymplectic manifolds) are not unique. 

In more detail, the manifolds Ti and Te can be covered with open patches such that each 
one on Ti is equivalent to one on Te, but the map cannot be extended to the outside of the 
patch, because it diverges at the boundary. Scattering trajectories contracting from right 
(left) and expanding to right (left) can both never lie within one and the same patch. Where 
two patches overlap, the two corresponding maps differ by a symmetry. The huge amount 
of symmetry in the Potsdam approach is the cause for the non-uniqueness of the map. 

Let us dwell a little more on these results. The local equivalence explains why the 
same radial equation results from both systems. What is, however, a possible physical 
significance of just local but not global equivalence of constraint dynamical systems remains 
unclear. On could imagine, for example, pasting together several copies of Ti and Te using 
the local equivalence maps. The result, however, seems to be a (possibly non-Hausdorff) 
presymplectic manifold with no reasonable physical interpretation (this is discussed at the 
end of Sec. |VTD . 

It is also conceivable that the global inequivalence of the systems is not important for 
the quantum theory; the delicate points where the map can diverge, lie all at the boundary 
between the scattering and the bound trajectories. One can speculate that this discontinuity 
is not relevant for the quantum mechanics, because the bound states become discrete anyway. 
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If we cut out this less dimensional boundary from both classical systems, we obtain systems 
that are globally equivalent. Further study is necessary 

Suppose next that some form of weakened equivalence between the two systems makes 
sense. Then, given an equivalence map, we can consider the Warsaw system as a part of 
the Potsdam one, and the missing information about the position of the shell in the right 
subspacetime can be provided. The Warsaw variables can be regarded as coordinates on 
a part of the Potsdam system. In fact, it turns out then that the Warsaw variables are 
coordinates of the shell in the left subspacetime only according to their name; in reality, 
they play the role of coordinates in the right subspacetime. Hence, the position of the shell 
in the right subspacetime is well-determined; paradoxically, it is the position of the shell in 
the left subspacetime that is uncertain (see the discussion at the beginning of Sec. |V|). 

Applying this point of view to the quantum mechanics of Refs. and f?| , we can make 
some progress. The self-adjoint extension of the Hamiltonian defines a unitary dynamics 
in the Warsaw coordinates of a part of the Potsdam system. It follows that each wave 
packet reemerges, during the expanding part of this unitary evolution, in the same right 
subspacetime in which it originally started to contract, because the variables describing it 
all the time are coordinates on this subspacetime. This is discussed in more detail at the 
beginning of Sec. [V|. 

On the other hand, the non-uniqueness of the equivalence map has an unpleasant conse- 
quence: although the self-adjoint extension of the dynamics seems to be unique, it is so only 
with respect to a particular choice of Warsaw variables. Different equivalence maps lead to 
different choices of the variables and these, in turn, to different dynamics of the Potsdam 
system. This is explained at the beginning of Sec. [V[ The difference is measurable, because 
the resulting scattering time delays are different. This suggests that one either has to look 
for some additional principle that could, together with the self-adjoint extension, lead to 
a unique set of time delays, or to look for some interpretation of the non-uniqueness (like, 
say, a loss of information) or to look for another way of dealing with the singularity. Future 
research may clarify the point. 
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II. DESCRIPTION OF THE SHELL DYNAMICS 



In this section, we briefly collect and round off some results scattered in literature making 
so the paper self-contained. 

A. The shell spacetime 

A sherically symmetric thin-shell spacetime solution of Einstein equations will be de- 
scribed in this subsection following closely Ref. ||. Consider two Kruskal spacetimes M.i 
and Ai 2 with Schwarzschild masses E\ and E 2 . Let Ex be a timelike hypersurface in A^i 
and E 2 be one in A4 2 . Let Ex divides A^i into two subspacetimes, M.\ + and M.\-, and 
similarly E 2 divides M. 2 into M.2+ and M. 2 _\ we chose fixed time and space orientation in 
the two-dimensional Kruskal spacetimes so that future and past as well as right and left 
are unambiguous; let then Ai 2 + and M.\+ be right with respect to hA 2 - and M.\-. Let Si 
and E 2 be isometric; then the spacetime A^x- can be pasted together with the spacetime 
M. 2 + along the boundaries Ex and E 2 . The result is a shell spacetime Ai s . As everything is 
spherically symmetric, only the two-dimensional Kruskal spacetimes are relevant. 

The observers are assumed to be in the asymptotic region of M. 2 +. Given a shell space- 
time, we shall often leave out the indices 1 and 2, having right (left) energy E + (.EL), shell 
trajectory E, and the right (left) subspacetime M.+ {M.-). Thus, A4 + = M. 2 fl M. s and 
M- = Mi n M s . 

The three-dimensional shell surface E carries the energy-momentum tensor T kl of the 
shell. This will be assumed in the form of ideal fluid 

T kl = (p + p)T k T l + jrf kl , 

where p is the surface mass density, p the negative surface tension, T k a unite timelike vector 
field (the three- velocity of the fluid) tangential to E, and 7^ is the metric induced on the 
shell from the surrounding spacetime. Let p = p(p) be the equation of state. 
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The spherical symmetry and the energy-momentum conservation lead to the matter 
equation 

A 1J + P + M = °' 

where A is the surface of the shell, A := 4ttR 2 . We choose one of the particular solutions of 
the matter equation, p(A), and define the so-called mass function M(R) by 

M(R) := A(R)p(A(R)). 

For example, the dust equation of state, p = 0, implies that M(R) = const, and each value 
of the constant defines a particular solution. 

The spacetime around the shell already satisfies the Einstein equations; thus, the only 
non-trivial equation still to be satisfied is the jump condition, the so-called Israel equation 
T2fl . It implies (for a derivation, see Ref. 0) the following two equations for the embedding 



functions T t (s) and R(s) of the shell in A4 e , where e is a sign, e = ±1 and s is the proper 
time along the radial generators of the surface. 

A) The radial equation: 

R 2 + V(R) = 0, (1) 

where 

vm- M2{R) E+ + E - ( E +- E -^ 2 i i to) 
V{R) - "Dp r WW ' ( } 

and 

B) the time equation (valid for all future oriented shell motions) 
where 

2E f 



F f := 1 



R 



S 



is the Schwarzschild function. 

Here, T± and R are the Schwarzschild coordinates of the shell in M.± (they are of course 
singular at the four horizons of the two spacetimes M.±) 

There are two types of solutions: bound and scattering; the scattering trajectories can be 
divided into expanding (R > 0) to the right and to the left, and contracting (R < 0) from the 
right and from the left (the usual notions of out- and in-going are not adequate to dynamics 
in the Kruskal spacetime). We stress that the four possibilities are unambiguous (see Ref. 
||): for each (scattering) shell spacetime, only one of these is realized in both subspacetimes 
Ai± simultaneously 

B. Warsaw approach 

In this subsection, we describe a short-cut approach (cf. e.g. J7| and [[J) to the shell 
dynamics. The first step of the approach is that only a subclass of the shell spacetimes is se- 
lected from all spherically symmetric shell spacetimes as described in the previous subsection, 
by fixing the value, E, say, of the Schwarzschild mass E- in their internal subspacetimes. 
Then the spacetime M.\ is a fixed, complete Kruskal spacetime and each shell spacetime 
defines a trajectory Si in it. This trajectory satisfies Eq. ([!]) and ([D with e = —1. In Ref. [f| 
a proof is given that each such trajectory belongs to a unique shell spacetime: the condition 
that the shell spacetime satisfies the full set of Einstein equations determines the mass E + 
of the right subspacetime Ai + and the trajectory R+(t), T + {t) in it up to a constant shift 
of T + (t). In this sense, the shell dynamics can be reduced to a dynamics of a particle-like 
object on a fixed two-dimensional spacetime M.\. In Ref. ||, the corresponding variational 
principle was shown to be uniquely determined (up to a coordinate-dependent factor in front 
of the super-Hamiltonian) if the super-Hamiltonian was required to be at most quadratic in 
momenta, and the value of the momentum conserved due to the time-shift symmetry to be 
the negative of the Schwarzschild mass of the external subspacetime Ai+, that is the total 
energy of the shell. 
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Let us describe this variational principle. We will leave out the indices ' — ' and '1' in 
refering to M._ and A4\ within the Warsaw approach. The coordinates x a , a = 0, 1 on 
M. play the role of the canonical coordinates of the system; the corresponding canonical 
momenta are denoted by p a , a = 0, 1. The action reads 

Si = J dt( Pa x a -Afd), (4) 

where t is an arbitrary parameter along dynamical trajectories, M is a Lagrange multiplier 
and the super-Hamiltonian C\ is given by 

d := l -F{g ab p aPb + M 2 ) - Wg ab p a £ b - l -W 2 . (5) 

Here, g ab (x) is the contravariant metric on A4 (observe that the supermetric Fg ab is flat), 
£, a (x) is the Schwarzschild time-shift Killing vector on A4, F = — <7 the potential W(x) 

is defined by 

wm ^ E .?m, (6 ) 

and M(R) is the mass function of the shell matter. 

Ref. M contains a proof that the variation principle (B) implies the radial equation ([[]) 
and the time equation (|3]) with e = — 1. The time equation (|3|) with e = +1 must be added 
by hand. 

We observe that the extended phase space is four-dimensional and that there is one 
constraint. Hence, the system has just one physical degree of freedom; this can be chosen 
to be, for example, the radius of the shell. 

As one easily verifies, 

Thus, £, a p a is conserved. The value of ^ a p a is —E + , where E + is the total energy of the shell, 
or alternatively, the Schwarzschild mass of the external subspacetime of the shell spacetime. 
Thus, p a cannot be homogeneous in the velocities x a ; this is the reason for the presence of 
the vector potential W£ a in the super-Hamiltonian C\. 
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To do calculations, we shall choose the Kruskal coordinates U and V for x a ; let us 
recapitulate some properties of these coordinates. The spacetime metric has the form 

ds i = J^EL dudv ^ 

K,e K 

where k = k(—UV), and the function k(x) is defined in the interval (—1, oo) by its inverse, 



K (x) := (x - l)e x ; 

this implies the identity 

k'(x) = . } , r . 

K{x)e K(x > 

The transformation to the Schwarzschild coordinates is given by 



R = 2Ek(-UV), T = 2E\n 

it follows that 



(7) 



m = -™ = !zl. (8) 
The time-shift Killing vector £ has the form 

i = ^ E {-ud u + vd v ). 

The variational principle @) and @ for the Warsaw approach reads in the Kruskal 
coordinates as follows: 

S x = J dt(PuU + P V V - AfCi), (9) 

where 

UVPuPy K-l W W> 

Cl ~ + ^T M + 4E {UPu ~ VPv) ~ (10) 

Our definitions (|^), ( [TDD an d @ °f ^i an d ^ differ from those given in Ref. by the 
factor F so that the super-Hamiltonian C\ and the potential W are smooth everywhere. 

11 



Using Eq. (|10|) one can show that the function C\ has a non-zero gradient everywhere on 
the phase space with the exception of the points with U — V — (whose projections 
to the configuration space is the crossing point of the two horizons R = 2E). Moreover, 
the constraint equation C\ — has no finite momenta solution there in the generic case 
W(2E) 0. This is roughly a consequence of the requirement that £ a p a = E + and the fact 
that all components of the vector £ vanish at the crossing point. Of course, in the special 
case of flat spacetime (E- = 0) this problem is avoided. We shall also show in Sec. |l[B| 
that this singularity is due to the choice of the coordinates x a and p a . 

An analysis that is based on the Kruskal coordinates cannot include the two special 
cases, E_ = and E + = 0, in which one of the subspacetimes is flat because the Kruskal 
coordinates become badly singular in the limit E± — > 0. Especially the case E_ = is, 
however, important. We give a short description of it in the Appendix. 



C. The Potsdam approach 

In this subsection, another approach to the shell dynamics is described. It starts so to 
speak from 'the first principles': in Ref. jRJ, general Lagrangian and Hamiltonian formalisms 
for massive shell are developed, starting from the Einstein-Hilbert and a fluid actions. The 
Hamiltonian formalism is reduced in Ref. JD] by spherical symmetry, using a transforma- 
tion of canonical variables invented by Kuchaf |13| . The basic properties of the resulting 
constrained system are as follows. 

The canonical coordinates involved are four coordinates of the shell in the left and the 
right subspacetimes Ai±. For example, one can take the Kruskal coordinates U± and V±. 



Another variable is the co-called 'Kruskal momentum' P% that was introduced in Ref. |LT 

dX e 



Plr '■= iLarctanh- 



K . ... ^ 

where X e := —U e + V e and Z t := U 6 + V e , so that the argument of the arctanh can be 
considered as the 'Kruskal velocity', that is the velocity of the shell with respect to the 
Kruskal frame (n e (l, l),n e (— 1, 1)) (n e is a suitable normalization factor). First, 
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observe that dV e /dU e > holds along each non-spacelike curve. There is a relation of the 
Kruskal momentum to the proper-time velocity (dU e /ds,dV e /ds); it follows from Eq. (|TT|) 
and the normalization condition, 

_(AE 1 fdU 1 dV 1 = _^ 
K e e Kt ds ds 

where K t := k(— U e V e ), namely, 



dU e _ v^e^ „-p k / R , 



e- F V H % (13) 



ds AE f 

* e e p y R <. (14) 



dV e _ y/K e e K < „p*/r, 



ds AE e 

The variational principle of the Potsdam approach has been written down in terms of Kruskal 
coordinates U e , V e , Kruskal momenta P^ and the energies E e in Ref. [11| : 

S 2 = Jdt {[-2EkP k + E 2 (k + l)e- K (VU - UV)\ + 2v[Ek] - vC 2 ) , (15) 

where 



Co : -- 



E^e- K ' 2 {-Ue PK ' R + Ve~ PK / R )} + M{E + k + + E„k_), (16) 



the symbol [/] denotes the jump of the quantity / across the shell, [/] := /+ — /_ and / 
denotes the average of the values of / from left and right: / := (/ + + /_)/2. 

Let us stress that the origin of the symplectic structure defined by the Liouville form in 
action fli~5|) lies in the form of the Einstein-Hilbert Lagrangian. As concerns the constraints, 
there are three of them in all: the two primaries, R + — i?_ = and C 2 = 0, and one 
secondary, \ = 0, which can be obtained by variation with respect to Pr± after employing 
the constraint E + k + = E_k_: 

dC 2 dC 2 



<9P+ dP K 
We obtain 
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X 



1 ( Ug Pk/R + Ve -W) 



The pair (x, R+ — R-) form the second-class part of the constraints set (see Ref. 

Hence, we have an eight- dimensional extended phase space and three constraints. This 
implies that there are two physical degrees of freedom; they can e.g. be chosen as £L and 
R_. This should be compared with the Warsaw approach of the previous subsection, where 
the Schwarzschild mass E_ was just a parameter having vanishing Poisson brackets with 
everything. 

The Kruskal metric is invariant with respect to the one-dimensional transformation group 
g x , A G (-00,00): 

U i-> Ue x , V h-> Ve~ x . (17) 



From this isometry, an infinite-dimensional abelian group of symmetry for the action ([151) 
can be constructed as follows. First, it can transform the coordinates LL, V_ of the left and 
U + , V + of the right subspacetimes independently. Second, this double transformation can 
be performed for each value of the pair (E + , EJ) independently. Such a transformation has 
the form 

U € ^ U e e x ^\ V e I— > V e e~ Xt ^ Ee \ e = ±1, (18) 

where X_(E_) and X + (E + ) are two arbitrary (smooth) functions. Let us denote the transfor- 
mation by G[X-(E-), X + (E + )]. Clearly, G preserves the products U £ V e and the two masses 
ER 



3 The origin of this symmetry is easy to understand: continuous symmetry groups are generated 
by perennials (constants of motion) ; i£_ and E + are perennials, and any functions A(E-) and 
A(E + ) are also perennials. Thus, one obtains a symmetry depending on two arbitrary functions 
of one variable. In fact, the symmetry group is even larger: it is generated by all perennials of the 
form A(E + ,E_)! 
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The corresponding transformation of the Kruskal momenta is determined by Eq. (|TT|): 

P € K ^P e K -RK(E e ). (19) 

Then, the functions U e e p i < / Re and V e e~ p K' Rt are invariants, and so are both constraints, C2 
and R + — A little more difficult is to see that the Liouville form of the action fll5"|) 
changes at most by a closed form. Let us write out the transformation of just the e-part of 
it leaving out the index e 

-RdP K + E 2 (l + K)e~ R (VdU - UdV) -> 

-RdP K + E 2 (l + K)e~ K (VdU - UdV) + Rd(RX) + 2E 2 (1 + n)e~ K UVdX. 

Using Eq. (^) and R = 2Ek, we obtain for the last two terms: 

Rd(RX) + 2E 2 (1 + K)e' K UVdX = RXdR + —dX + 2E 2 —dE, 

2 dE 

which is a closed form. 

The Kruskal spacetime also admits an isometry group of four elements with generators 
T\ and T2, the two transformations that are obtained by extending the Schwarzschild time 
reversal T 1— ► — T, either from the quadrant U < 0, V > or from U > 0, V > to the 
whole of the Kruskal manifold: 

71(17, V) = (-V, -U), T 2 (U, V) = (V, U). 

These transformations can be used to define symmetries of the variational principle (|H|) as 
follows. Consider first 7\; it changes the time orientation of the Kruskal spacetime. Thus, it 
must act in both Kruskal subspacetimes of each shell spacetime simultaneously, or else the 
resulting shell spacetime has no time orientation. We have, therefore 

Ti{U e ,V e ) = (-V e ,-U e ). 

The action of this transformation on the momentum P^ is determined by Eq. (|TT|): 

Pk - -Pfc- 
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Finally, E £ are clearly invariants. This defines a transformation in the extended phase space 
that will be also denoted by T\. Observe that the constraints C2 and R + — R_ are invariant, 
but that the Liouvile form in the action (|T5| ) changes sign. Thus, T% is an anti-symplectic 
map (this is common for time reversals). 

Next, consider T 2 ; it changes the space orientation exchanging right and left. Thus, we 
have to define: 

T 2 (U e ,V e ) = (F_ e ,l7_ e ), 

and, similarly, 

One can see immediately that both the constraints and the Liouville form are invariant with 
respect to this transfromation. 

III. CARTAN FORMS 

A. Reparametrization-invariant reduction method 

The two variational principles (^) and (|T5D seem to be as different as one can imagine. 
Although they depend on some common variables (because U = U- and V = V_), the 
momenta are different in each case, and the nature of the super-Hamiltonians fllCf ) and (|16D 
is very different: C\ = is a condition on the norm of the two- velocity of the shell in the left 
subspacetime, whereas C2 = relates the time-time-components of the second fundamental 
forms of the shell in the left and right subspacetime (indeed, C2 = is a component of the 
Israel equation). 

It turns out that the simplest way to compare the two systems is to solve all constraints 
thus reducing the systems completely. It is, however, the kind of reduction that does not 
deparametrize the system: no choice of gauge is necessary. The resulting Lagrangian is the 
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value Q(X) of the so-called Cartan form G at the tangent vector X of a trajectory, and so 
it is homogeneous in velocities. The variation of the corresponding action, 

s = J dte(x), 

leads to dynamical equations of the form 

d&(X,5X) = 0, V5X. (20) 

dQ coincides with the presymplectic form •) on the constraint surface of the system, and 
the equation of motion ( ^0|) simply says that X must lie in the degeneracy subspace of Q 
(for more about presymplectic forms see, e.g. Refs. [14]] and [15], cf. also ||16|| ). 



B. Warsaw description 

In this subsection, we will reduce the Warsaw system. Since the Potsdam approach was 
written down only for massive shells as yet, we restrict ourselves to massive shells, M(R) ^ 0, 
in the rest of the paper. We observe then that the constraint C± = can be easily solved, if 
we express the momenta P\j and Py by means of Pk as follows. The variation of the action 
(fj) with respect to the momenta gives their relation to the velocities: 

2(2£) 2 (U WU\ 



P 



u 



UV \N AE ) 
2{2E) 2 (V WV\ 



P V = -±7+ -T7 + 



UV \J\f 4E J 



If we insert these expressions for the momenta into the constraint (]1Q|) , we obtain that 

U V _ (k - 1)V M 2 
Jf'Jf~ k {4E) 2 ' 

A comparison with Eq. ( |T2"D and a use of Eqs. (|T^) and ( |H] ) leads to 



f M P m W\ 
P u = 2E{-^=e p «/ R + -), (21) 



y/Ke* U 

JL e -P K /R _ ? 
^/k^ V 
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One easily verifies that these expressions for Py and Py satisfy the constraint (|T0|) identically. 
In this way, we have arrived at the constraint manifold I\ with the coordinates U, V and Pr. 
On this manifold, there are two important quantities: the pull-back —E + of the conserved 
quantity £ a p a = 1 / '(4£7)(— £/ Py + V Py) and the Cartan form. Substituting for the momenta 
Py and Py, we obtain: 

E ^ E + ^B { - UeP '" H + Va ""^ ~ ^HT' (23) 

where R is an abbreviation for 2Ek(—UV). The Cartan form is obtained, if we substitute 



for the momenta into the Liouville form of the action 
0i = 2EWd fin 



2E ^ffl U p «l R dU + e~ p «/ R dV) . (24) 
\/Ke K v 1 



One would expect that these objects would be singular at the point U = V = where 
the horizons cross, and indeed, the Cartan form (|24] ) diverges even at both horizons U = 
and V — 0. However, this singularity can be removed by subtracting the differential of the 
function 



2EW{2E) In 



in this way, we obtain an everywhere regular Cartan form 

lr = 2EW r (VdU - UdV) - ^£ (e Pli/R dU + e~ Pl</R dV 

where 

m _ uv) ;= W(R)-W(2 E) 

is a smooth function. For example, with M = const, 

M 2 



W r 



AEne K 



The regularity of the constraint system requires, however, a slightly stronger condition 
that just the smoothness of the Cartan form: the presymplectic form dQ± = d®\ r must 



be smooth and its degeneracy subspace must be one- dimensional everywhere in T\. Let us 
calculate the presymplectic form; a straightforward but tedious procedure yields: 

+ MP « ( Ug Pk/R + Ve -PK/R) )dUAdV 
K 3 e K V ne K v ' J 

] ' (e Pli/R dP K AdU- e~ PK/R dP K A dV 



where the indices R and k denote partial derivatives with respect to the corresponding 
variables. An inspection shows that this form is smooth and non-zero everywhere; this is 
sufficient for a two form on a three-dimensional manifold to define a smooth one-dimensional 
degeneracy distribution. Thus, the constraint system (r 1; Qi) is completely regular even at 
U = V = 0. The singularity of the variational principle @ is due just to the choice of the 
variables x a and p a : the simplicity of the super-Hamiltonian (|5]) is traded for the singularity 
at U = V = 0. In any case, the constraint manifold T\ is defined by the coordinates U, V, 
and Pk in the ranges 

-UVe(-l,oo), £+e(0,oo). 



C. Symmetry and adapted coordinates 

In this section, we shall study the static symmetry of the Kruskal spacetime and the 
corresponding symmetry of the system (Ti,Oi). We shall use the symmetry to find some 
coordinate systems which will simplify the subsequent calculations. 

Eqs. fll7), and ( |19D for A = const, define a transformation on I\; this transformation 



leaves the function E + and the form 0i invariant. Thus, it is a symmetry of the constraint 
system. From the discrete transformations, only T\ survives; T2 would map our system to 
an equivalent one, which would be obtained if we based the description on the trajectory 
in the right subspacetime M.2- It is quite natural that this much smaller transformation 
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group remained from the group G[\+(E+), A_(£L)] here, because the right subspacetime is 
missing and E = const, so X(E) = const. We call this transformation G(\). 

The orbits of the group G(\) in I\ are one-dimensional, and the quotient Ti/G(X) is 
a two-dimensional manifold. It is useful to introduce coordinate systems that are adapted 
to this quotient structure: two coordinates that are constant along the group orbits, and 
another one that transforms simply by the group. 



1. Coordinates u, v and P 

Let us define: 



Ue P K /R Ve -P K /R 
U-= ,==- , v:= — =-. (25) 



The factor 1 / \J 'ne K does not improve the transformation properties, but it leads to additional 
simplifications: indeed, all exponentials, square roots and transcendental Kruskal functions 
will disappear. Futher, 

P := (26) 

hence 

P h-y P- EX. 



Eq. (Rl) implies that 



or 



uv = -l + 2 -f, (27) 



2E . . 

R = ——, 28 
1 + uv 



and this, together with Eq. (0) yield 

1 



(29) 

1 + uv v ' 
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The objects E + and 0i can be transformed into these coordinates with the results: 

„ M(R). . M 2 (R) , . 

E + = E- ~ (u -v)- (30) 



and 

F,M 2 (m / it \ 

-2EM(R)d(u + v) 



* = ( in 



M(R) <u +v)(l+ ^ ) d/2 - 4£+dP, (31) 



2 y ' \ R 
where R is defined by Eq. (p8|). 



Coordinates R, E + and P 

The function I£ + is not only invariant with respect to the group G(X), but also an integral 
of motion. Hence, projections of the dynamical trajectories to the quotient Ti/G(X) are just 
the curves E + = const. It will prove advantageous, therefore, to transform to the coordinates 
R and E + on the quotient. Of course, R and E + are not everywhere regular coordinates: 
clearly, at the maximal radii of bound trajectories, gradients of both functions R and E + 
vanish in the direction of the trajectories; we can, however, still work with these coordinates 
in the rest of the quotient, and match the results at the singular points. The transformation 
from u and v to E + and R are given by Eqs. (|30|) and (28) . The inverse transformation is 
obtained by solving these equations for u and v: 



-Y + luVY 2 - AX 
Uu ~ 2 : 

Y + tuVY 2 - AX 



where u is a sign, whose significance will be established shortly, 

\ M(R) 2R J ' 

and 
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(32) 
(33) 



Y 2 - AX = -AV(R), (34) 

where V(R) is given by Eq. @. Thus, if V(R) vanishes, R must also vanish and this happens 
only at the turning point of a bound trajectory. The meaning of the sign u in front of the 
square root is simple. As 

(-it + v) u = Y, 



(u + v) u = uVY 2 - AX, 

it distinguishes the upper and lower half of the w-f-plane; the turning points are lying at 
u + v = and the coordinates E + and R are not regular there. 

D. Potsdam description 

In this subsection, we reduce our second system completely. In analogy with Eq. 
we define 

UeB PyRe V e e- p y R ° 
u £ := . = , v e := — , (35) 

and obtain the relations 

2E 6 2E £ 1 

u e v e = -1 + — , R € = — , K e = — . (36) 

R e 1 + u e v e 1 + u e v e 

The three constraints Ci = 0, R + = R_ and x = simplify greatly in these coordinates: 

C 2 = ^[R(-u + v)]+M(R), 

! r i 
X = 7,[u + v\, 

R = = R-. 

They can be solved immediately, either for u + , v + and E + : 

u + = u_ + «g>, (37, 

» + --^, (38) 

M(i2) . , M 2 (i?) 

E + = E„- —^(u- - «_) - (39) 
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or for it_, V- and P_: 



P+ + — ^-(u- 



M 2 (P) 



2/? 



(40) 



(solutions for U- and u_ are given by Eqs. fl3~T|) and (|38|)). 
From the definitions Eqs. (B5|), we obtain 



In 



Vr 



IP. 



K 



R 



+ In 



(41) 



To calculate the Cartan form 02, we first express the Liouville form of the action (|15D 
in the variables w_, u_, P^, P_, u + , t> + , P^ and P + using an analogy of Eq. (§) 

M e w e = (1 - K e )e Kf ; 

after substituting for E + , u + and t> + and employing several times the equations R = 
2P_k_ = 2E+K+ we obtain: 



e 2 = - 



^--E 2 )d[\n 



+ 4EdP 



(42) 



where the sign '=' suggests that we have added some closed forms, and 

E f Pl 



P :-- 



J ^ K 



R 



(43) 



The form is well-defined and smooth everywhere on the constraint manifold T 2 that is defined 
by the ranges of the coordinates it_, t>_, E_, P_ and P + : 

■u_t>_ G (— l,oo), u + v + G (— l,oo), P_ G (0, oo), 

E + e (0, oo), P_G (-00,00), P+G (-00,00). 
The reduced action has the form 



S2C — I dtL 2 c, 



where the Lagrangian is 
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L>2C — — 



R 2 



— E 



u v 



(44) 



The symmetry group of the system (T 2 , G 2 ) contains the whole infinite-dimensional group 
G[A_ (.EL), X + (E + )]. In the coordinates E_, w_, t>_, P_ and P + , this transformations take 
on the from: 

P_^P_, v-^V-, P_ ^ P_ - ELA(EL), P+ ^ P+ - E+\(E+). 

2 changes by AE + d(E + \ + (E + )) — 4P_d(P_A(P„)), which is obviously a closed form. From 
the discrete group, 7\ remains simple; it acts as follows: 



P_ I— > E_, U— I — ^ — M_, f_ I — ^ — P_ 1-^ — P_, P + 1-^ — P 



+■ 



IV. EXCLUSION OF A CYCLIC COORDINATE 

If one compares the systems (T^Oi) and (r 2 ,6 2 ) the first difference that catches one's 
eye is that T 2 is five- and I\ only three-dimensional. In this section, we shall get rid of two 
dimensions of T 2 . The fact that the variable E appears in (I^, Q x ) as a parameter implies 
that we have, as the first step, to select the submanifold T 2 e C T 2 of constant coordinate P_ 
in T 2 . As E_ is a constant of motion, the dynamical trajectories that intersect T 2 e all remain 
in T 2 e- However, the dimension of Y 2E is still larger by 1 than that of T 2 ; moreover, the 
pull-back Q 2E of 6 2 to Y 2E defines a presymplectic form dQ 2E that has a two-dimensional 
space of degeneracy: it contains the vector <9/<9P_ in addition to the direction of motion. 
Thus, in the second step, we have to construct the manifold T E := T 2E /P_; it is the quotient 
of T 2E by the P_-curves (i.e. curves with w_ = const, t>_ = const and P + = const). A form 
Q' 2 e that differs from Q 2E by a form closed on T 2E does not contain the variable P_; it 
can, therefore, be pushed forward by the quotient projection ir to T E giving a form that we 
call Qe- d@E is a presymplectic form on with a one-dimensional degeneracy subspace 
everywhere, and the integral manifolds of this subspaces coincide with the projection to 
T E of the original dynamical trajectories. Observe that we have to take a quotient so that 
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P_ remains arbitrary; we are not allowed to choose a surface transversal to the P_-curves 
instead, defined for example by the equation P_ = const, for P_ is not a gauge coordinate 
or an integral of motion (it satisfies a non-trivial dynamical equation that would be violated 
by P_ = const). 

In terms of coordinates and dynamical equations, the construction is very simple. The 
pull-back Q 2E of 6 2 to Y 2E coincides with 9 2 in the coordinates w_, i>_, P_ and P + ; only 
P_ changes its status: it becomes a constant parameter. Then the term 4P_<iP_ is the 
differential of the function 4P_P_ on T 2 e and can be omitted because it does not contribute 
to d0 2 E'i thus, we end up with the form 



this form is independent of P_ and so can be pushed forward to Te] in the coordinates w_, 
f_ and P + , the push-forward is given by the same expression as Q 2 e- Thus, the new action 
reads 



on the space T E with coordinates w_, t>_ and P + . Let us denote this three-dimensional 
system by (Te,@e)- The symmetry of this system is still infinitely-dimensional: it is the 
transformation group G[\+(E + ), A) (the e = +l-part of the transformation acts only on P + 
and A_(P_) is a constant, because E_ is) and the time reversal T\. 

Let us compare the dynamical equations obtained by varying the actions S 2 c and Se- 
The following observation is helpful: the part of the Lagrangian ([0|) that contains the 
variables w_, t>_ and P + coincides with the Lagrangian (pp. It follows that the variations 
of action S 2 c with respect to the variables w_, t>_ and P + differ from the corresponding 
variations of the action Se only by terms that are proportional to P_; but P_ = is the 




(45) 



dtL 



where the Lagrangian is 




(46) 
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dynamical equation obtained by varying the action S2C with respect to P_. Hence, the 
three dynamical equations of the system (r 2 , 2 ) due to variation of w_, t>_ and P_, if one 
sets P_ = in them, are identical with the corresponding three equations of the system 
(T E ,0 E ). The equation E_ = is free in the system (T E ,Q E ), because E_ is a constant 
parameter there. Finally, there is another equation in the system (r 2 , 2 ), namely that due 
to variation of E_; it can be written in the form 

4>- = (47, 

it cannot be obtained in the system (T E , Q E ), because P_ is the cyclic coordinate that has 
been eliminated, and the new variation principle does not contain any information about it. 
We observe that P_ is 'completely smeared' even classicaly, as any point of T E is a whole 
P_-curve, P_ G (—00, +00), and the remaining variables in the system (T E ,Q E ), namely 
U-, V- and P + have vanishing Poisson brackets with P_ in the original Potsdam system. 

This is in a nice correspondence with the quantum mechanics construted for the two 
classical systems (if the factor-ordering problem is solved suitably): the Hilbert space TC of 
the system (r 2 ,02) can be written as the orthogonal sum of the eigenspaces Ti E _ of the 
operator P_; P_ is a point of its spectrum cr(P_). On the Hilbert space 7i E _, only those 
elements of the algebra of observables can act that commute with E-, and TL E _ is the Hilbert 
space of the system (T E ,Q E ). 

We can even reconstruct a trajectory of the system (r 2 , 2 ) from one of (r^, 9^) using 
the equation ( $7\ ) which must simply be added by hand, as follows. Let V-(t) and 

P+{t) be a trajectory of (r^, 0#) parametrized by an arbitrary parameter t. Let us substitute 
these functions for w_, t>_ and P + in Eq. ([47]). We obtain an equation of the form 

dP— 1 dL Ej , , , , . - . ,. 

the integration will yield the function P_ (t) depending on one arbitrary constant (this con- 
stant is determined by an initial value of P_). Such a procedure is, in fact, equivalent to 
reconstructing the manifold T 2 from by 
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r 2 = T E X (-00,+00)£;_ X (-00,+00)p_ 

and adding the term AE_dP- to Qe- In the quantum mechanics, this corresponds to defining 

(in general, one has to be careful about the spectrum of the conserved quantity; this is here 
just (— oo,+oo)) and 

~ i d 



the independence of the Hamiltonian of the system (r^, Qe) will then lead to some time 
evolution of P_. 

V. SEARCH FOR THE EQUIVALENCE MAP 

It remains to find a transformation between the system (T El Q E ) and (r^Gx) showing 
that they are equivalent. This problem will be studied in the present section. 

The basic properties of such a map (which is, in fact a 'morphism of presymplectic 
manifolds'), let us denote it by $, are the following 

1. $:rii— ^r^isa diffeomorphism. 

2. The pull-back $*0£ of Qe to I\ differs from 0i by a closed form. This guarantees 
that the presymplectic forms coincide^. 

3. If fx : Y\ i— ► R and /e : k> R are quantities with the same physical or geometrical 
significance in both systems, then = fx- 



4 To define equivalent dynamics, two presymplectic forms just had to be proportional to each 
other with an arbitrary function on the constraint surface as the factor of proportionality. However, 
presymplectic forms contain more information: they define the Poisson brackets of perennials (that 
is, functions that are constant along the trajectories, cf. |1~6[| ); thus, the factor must be equal to 1. 
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Suppose that such a map exists. Let us express it by means of the coordinates (u, v, P) 
on (Ti,0i) and (u-,v_, P + ) on (Te,®e) '■ 

u_ = u_ o v,P), (48) 
v _ = v _ o $(u,v, P), (49) 
P + = P + o$(u,v,P). (50) 

The functions ©-© can be viewed as a coordinate transformation and u, v and P as 
new coordinates on Te- This is an interpretation of $ that avoids the following paradox: 
P_ is to be totally undetermined on Te (see the previous section) and at the same time 
P_ has the same physical meaning as P, which is to be a well-defined function on Tg by 
Eqs. fl4"8"|) -([50|) . By considering Eq. (|50|) as a coordinate transformation, we regard P as a 
coordinate on the right subspacetime, because P + and P are in one-to-one relation [see also 
Eq. ( [52"D and fl5"3|) l and P + is such a coordinate. 



Another important point is that there is an infinity of different maps that satisfy above 
requirements, if one exists. This follows from symmetry: given a map $, we can sandwich it 
between the symmetries to obtain another one. The non- uniqueness of $, however, leads to 
an uncertainty in the coordinate P + (and so to that of the Schwarzschild coordinate T + ) of 
the shell. The symmetry G is, from the point of view of the spacetime geometry of each shell 
spacetime, nothing but a remainder of the original general covariance: the Schwarzschild 
or the Kruskal coordinates are not uniquely determined by the geometry of a fixed Kruskal 
spacetime. To map the dynamics of the system (ri,Oi) to (Te,@e), a particular $ must 
be fixed. Which one does not seem to matter as long as (classical) physical properties are 
concerned: a change in $ amounts just to relabeling the classical dynamical trajectories. 

In a quantum theory such a relabeling may lead to a problem, however. As an example, 
consider the unitary extension of the quantum shell dynamics described in Ref. 0. There, 
one works with coordinates analogous to u, v and P and one uses the simplicity of the 
Hamiltonian in these coordinates to show the existence (if the rest mass of the shell is 
smaller than Planck mass) and uniqueness of a self-adjoint extension of the Hamiltonian. 
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The eigenfunctions of the extended Hamiltonian are linear combinations of contracting and 
expanding waves. One can construct wave packets from them that are, at each time T, 
spacially concentrated around a well-defined wave-function maximum so that one can plot 
the radius Rm of the maximum as a function of the time T. The function Rm(T) diverges 
for T — > — oo, then decreases till, at some T , a minimum R M (T ) > of R M {T) is reached, 
and then increases again to infinity as T — > oo. For each wave packet one can define a time 
delay AT between the departure of the packet at R = oo and the arrival at R = oo. Of 
course, AT is not the limit \im R=00 [T 2 (R) — Ti(i?)], where Ti i2 (R) are the two solutions of the 
equation Rm(T) = R, because the difference T 2 (R) —Ti(R) diverges as R — > oo; one defines 
AT e.g. by comparing the packet 'trajectory' with some standard scattering trajectory. 

To interpret this scattering, we need, however, the function R M (T + ) rather than R M (T). 
In principle, we can calculate the first from the second by using some fixed map $. The prob- 
lem is that two different maps, $1 and $ 2 , say, will then lead to two different Rm{T+) that 
imply in turn two different time delays AT. This is clear because contracting (expanding) 
part of Rm(T) lies in the subset of the phase space that contains contracting (expanding) 
scattering trajectories and the difference between $ x and $ 2 along the expanding trajectories 
is independent from their difference along contracting trajectories. 

A. The differential equation for $ in the general case 

In this subsection, we shall return to the classical theory and reformulate the above 
requirements on $ in the form of a differential equation. We assume that E > and 
E + > 0; the special case E = is studied in the Appendix. We have two manifolds, r\ with 
coordinates u, v, and P that satisfy the conditions R > and E + > 0, and T E with the 
coordinates U-, i>_, P + satisfying the same condition (which is independent of P_ or P + ). 
According to the point 2, 

d($,e E -e 1 ) = o. (5i) 
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The functions u and w_ as well as v and t>_ represent the same physical quantities, so 
according to the point 3 the transformation (f48|)-(|50|) has to preserve them: 



it_ o &(u, v,P) = u, (52) 
v_ o $(u, v,P) = v. (53) 

Thus, the only non-trivial part of $ is the transformation (ptl; let us denote the correspond- 
ing function by </>: 

<j>(u,v,P):=P+o$(u,v,P). (54) 
Substituting Eqs. (0)-® into Eq. (|5T]), we obtain 

d (®E\ E _ =E>u _ =u , v _ =v , P+=<t> - 6x) = 0. (55) 

This is a partial differential equation of the first order for the function 0; its characteristics 
are easily found to coincide with the curves E + = const. We can, therefore, reduce this 
equation to an ordinary differential equation, if we transform it to the coordinates R and 
E + . That is our next task. 

Let us define the functions A, B, A' and B' of R and E + by 

®e\ e _ =E!U _ =u , v _ =v , P+=4 , = AdR + BdE + - 4E + d<P, (56) 
6! = A'dR + B'dE + - AE + dP. (57) 



Substituting Eqs. ( |56|) and (|57D into Eq. ([55"D, we obtain the equations: 

Ae + ~ A' E+ - 4<t> R = B R - B' R , 

0P = 1 

(the indices denote partial derivatives). These two equations are equivalent to the system 

= P + A(R,E + ), (58) 
1 

dR 4 ( 



^ = \(A E+ -B R -A' E+ +B' R ). (59) 
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Eq. ([59]) is the desired ordinary differential equation (for the function A). 

Let us work out the explicit .R- independence of the right-hand side. We obtain from 



Eq. and (|45|): 

A = 

B = 

A' = 
B' = 

We find easily 



- E 2 ) ( In 

- E 2 ) ( In 



r [ 4 



)(- 


u + 







E+ 



EM 2 

R 
EM 2 

R 



- 2EM(u + v) R -^-(u + v) (l + ^p) , 



E+ 



2EM{u + v) E+ . 



A E+ ~B R = 2E+ ( In 



u 

V 



In 



U- 



In 



M ( 2E 

- T (l + -](.M -V 



+ 2EMr(u + 



Next, we have to express everything in terms of R and E + , using Eqs. 
first step, we derive the following helpful equations: 

E + — E 

u = - 



and 



E 4 



M 
- E 



M VP 

h w , 

2i? Mi?' 



M 



M VV 
+ 2R + UJ MR> 



V 



-E M 

1 Yuj- 

M 2R MR' 



E + -E 
M 



M 

h LO 

2R 



V 

MR' 



where 

V(R) := -M 2 (R)R 2 V(R). 
A straightforward calculation then gives 



). As the 
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A E+ -B R =-^=(R + 2E + ) ( E+ -E + ^Pj + 

I 1 I L'- It/ L~ \ I ) O I / - L j 



M 



((4£^ - 4E + E)R + 2M 2 E-i 



A' - B' B = ^— f 2(£+ - -E)i? 2 + (4£ + £ - AE 2 + M 2 )R - 2EM' 

+ 2RVV V 



Collecting all pieces, we arrive at the differential equation for A in the form 
' /A f -2(E + -E) 2 + M 2 (E+ + E) 1 



dR 4 VP V v ^ 7 v ^ ' i? 



+ ((£+ " + M 2 (E + + E)J . (60) 



Observe that the upper and lower -u-f-plane give the same curves, only their orientation is 
different. 



B. Solution of the equation for A in the case of dust 

In the general case, we cannot say much about the solution of Eq. (J@Q), because the 
function M(R) in it is an arbitrary (positive) function. In the case of dust, M = const, 
however, the equation is readily solvable. In this case, we obtain 

% = M- 2{E +- E? + m2(E + + E)1 b)- (61) 

Eq. (|6l|) can be solved by elementary integration; the solution, however, will change its form 
if the parameters E + , E and M vary. Thus, we obtain only a local form of the function A. 
In the present subsection, a careful discussion of the different cases is given. The results of 
this discussion will be used in the next subsection where we try to match the different cases 
smoothly together. 

Let us write the solution of Eq. (^T[) in each of the half u-v-pl&ne, u = ±1, as follows: 

A w = A£ + A^ + A£, 

where Aq = Aq(E + ) is an arbitrary function of E + (integration constant — this is the non- 
uniqueness in $ due to the symmetry), 
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au . , . (E + -Ef r dR 
-~ U 2 J 



and 

. M 2 (E + + E) /■ dR 

For dust, the function V becomes a quadratic polynomial of R: 

V(R) := ((E+ - Ef - M 2 ) R 2 + (E+ + E)M 2 R + 

The discriminant of the quadratic equation V(R) = equals to M A (AEE + + M 2 ); it is 
always positive, so there are two roots: 

M 2 1 
2 (E+ + E) + yJAE+E + M 2 

and 



M 2 (E+ + E) + yAE + E + M 2 
R2 = ~~2 (E + - Ef - M 2 • (62) 

R 2 is positive for (E + - E) 2 — M 2 < 0, that is E + G (Max(0, E - M),E + M). Then 
V(R) > in the interval R G (0, R2) and the trajectory is bound, with maximal radius R2. 
For {E + - E) 2 - M 2 > 0, R 2 < i?i < 0, > in the interval i? G (0, 00) and the 

trajectory is unbound. 

To calculate the function A^, we have to distinguish three cases. 
(E + — E) 2 — M 2 > 0. In terms of the variables E, E + and M, this means that either 
E + > E + M or E + < E — M. The corresponding trajectories are the contracting and 
expanding scattering states. Then 

\ lj / n n \ (E + -E) 2 2aR + M 2 (E + + E) 

Af(R, E + ) = -u> y + _ ' arccosh v + 4, 63 

1 2^ M 2 JAE + E + M 2 



where ct := (E + — E) 2 — M 2 ; we use the increasing branch of arccosh. A^ is regular for all 
R G (0, 00). The values at the end-points are: 

A? (0, E + ) = - J E+ -^ r arccosh E+ + E 



2 V® ' ^ Jae+e + M 2 ' 
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and Af diverges logarithmically with R — > oo to — uoo. 

(E+ - Ef — M 2 = . Then i?i < 0, R 2 -»• oo and the polynomial > for all i? G 

(0, oo). The trajectories are the parabolic scattering states separating the scattering states 
from the bound states. The integral is: 



(E+-E 



,2 



Ar = 2(E + + E) ^ A{E+ + E)R + (64) 

It is again regular for all R G (0, oo). The values at the end-points are: 

*„, ^ x (E+-E) 2 M 
1V +; 2(E + + E) 

and it diverges as R}l 2 for R — > oo. 

— £) 2 — M 2 < 0. In terms of the variables E, E + and M, this means that E + G 
(Max(0, E — M), E + M). Then R 1 < < i? 2 - The corresponding trajectories are the bound 
states and the integral is 

(E+-E) 2 2aR + M 2 (E + + E) , . 

A?(R, E + ) = -uj k + - ; - arccos v + - 65 

2^ M 2 ^AE + E + M 2 

(we use the decreasing branch of arccos). Af is regular for all R G (0, R 2 ). The values at 
the endpoints are 

A^0,E + ) = - J E+ Zl )2 arccos , ^ + + E , 
2 V-a y^E+E + M 2 



and 



A?(ifc,^ = -o; ( ^ + J^ tt. (66) 
The function A 2 can be calculated immediately with the result 

2 2R X AE + E + M 2 



which is valid for all values of -E 1 , M, and R for which "P > 0. A£ diverges at R — > for 
all trajectories. For scattering trajectories, 

A, oo,£+ = — c<j ln — 

2 2 ^E~E^W 2 
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is finite. For bound trajectories, 



E+ + E (E+ + E) + JAE+E + M 2 
A%(R 2 , E+) = -u + In — ¥ (68) 



is also finite and for the limit R2 — > 00 (E + E ± M), we obtain 

IE i M. 

lim A%(R 2 ,E + ) = -uj ln(2£±M). 

R2=oo 2 



C. Matching and patching 

The function A(R, E + ) must be well-defined and smooth for all values E + > and R > 
of its variables, and this at all values E > and M > of its parameters in order that it 
determines a map $ with the required properties. The 'pieces' of A obtained in the previous 
section must, therefore, be smoothly matched together. In the present subsection we prove 
that this is impossible and give some discussion of the negative result. 

Within the half planes u+v < (oj < 0) and u+v > {uj > 0), the function A% (R, E + ) is 
smooth everywhere, but A( as a function of E + seems to be divergent at (E + -E) 2 -M 2 = 0. 
A routine inspection in the complex plane reveals that is, in fact, a smooth function 
there for each u. Hence, the sum Af(R, E + ) + A^(R, E + ) is smooth inside each half plane. 
The whole function A w can, therefore, be made smooth by an arbitrary smooth choice of 

^0 • 

The main problem is the matching at the boundary u + v = between the two half 
planes. Let us first study some properties of the curves E + = const in the w-f-plane. 
The Kruskal coordinates U and V are both future oriented; hence the past singularity 
R = lies in the quadrant U < 0, V < 0, and the future one lies in the quadrant U > 0, 
V > 0. The transformation ( p5[) [or (|35)] preserves signs, so the past singularity lies in the 



lower, and the future one in the upper half of the w-f-plane. Thus, in the lower half plane 
(u + v < 0), R is increasing along all trajectories and in the upper half plane (u + v > 0), 
R is decreasing. Observe that only bound trajectories can cross over from one half plane to 
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the other; scattering trajectories are always imprisoned inside one half plane: the expanding 
(R > 0) in the lower and the contracting (R < 0) in the upper half plane. 

Let us try to extend A continuously across the boundary u + v = using the remaining 
freedom in A%(E+). We shall split A%(E + ) into two terms, A% 1 (E + )+A% 2 (E + ). In the lower 
half plane, A{~ is increasing with R if E + is kept constant. If E + G (Max(0, E — M), E + M), 



Af assumes the value (p6 ) for oo = — 1 at the boundary u + v = 0. In the upper half plane 
Aj" is decreasing with R and it reaches the boundary with the value (|66|) for w — +1. Thus, 
we can construct a continuous function by choosing A§ X {E+) = and 

A+( B+ ) = ^±i££ I r. 



The function Af(R, E + ) +Aq 1 (E + ) in the upper half plane, together with A^(R, E + ) in the 
lower, define in fact a smooth function at all points the lower half plane, at the boundary 
u + v = 0, and in the subset E + G (Max(0, E — M), E + M) of the upper half plane; let 
us denote this set by X>_. T?_ consists of all points at the expanding scattering and at the 
bound trajectories. However, there is no continuous extension of this function to the rest 
of the upper half plane, because it diverges at the points of the upper half plane satisfying 
{E + — E) 2 — M 2 = 0, that is, at the contracting parabolic trajectories. 

Let us turn to the function A^(R, E + ) starting again in the lower half plane. It is 
increasing from the value — oo at R = to a finite value at R = oo along all scattering 
trajectories, and to a finite value ( |68| ) with uj = — 1 along bound trajectories at the boundary 
u + v = 0. The function QB8D has finite limits as E + approaches the values E ± M. Thus, 
there is no problem to extend this A^iR, E + ) to the whole upper half plane. One has just 
to choose Aq 2 (_E' + ) = and 



E + + E + JAE+E + M 2 
A+(E + ) = (E + + E) In -t ^ 

everywhere. The two functions A 2 (R, E + ) in the lower and A 2 (R, E + ) + Aq 2 {E + ) in the 
upper half plane form together a smooth function in the whole u-f-plane. 

It is also clear that no continuous choice of Aq X + Aq 2 can remove the singularity in the 
upper half plane. Hence, any allowed choice of A that is continuous in the lower half plane 
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will necessarily diverge at all points of the contracting parabolic trajectories (E + = E ± M) 
in the upper half plane. This result already means that there is no function A that would 
be continuous everywhere. 

An analogous construction starting from the upper half plane leads to analogous results: 
A i can be made smooth only in the open subset T> + of the w-f-plane that contains points of 
all contracting scattering and all bound trajectories. It diverges at all expanding parabolic 
trajectories. A2 can again be chosen smooth everywhere. Of course, the two constructions 
end up with two different solutions: in D_nD + , they differ by a function of E + that diverges 
at all points satisfying E + — E ± M. 

The two constructions in the foregoing paragraphs deliver, in fact, two maximal con- 
tinuous extensions of our local solutions for A. This two maximal continuous extensions 
are of course also smooth (if Aq are chosen so), but it is the continuity that is lost at the 
boundary. Similar conclusions can be drawn for the case E = or E + = 0, but the proof 
must use different vartiables, because u and v can serve as two regular coordinates only for 
E > 0. The proof for the special case E = is sketched in the Appendix. 

As discussed at the beginning of the present section, the map $ (if it existed) can be 
considered as equivalence map (morphism) between the two systems (Ti, ©1) and (T E ,Q E ) 
or, alternatively, as a definition of new coordinates, u, v, and P, for the system (Te,®e)- 
The fact that $ exists only locally and is non-unique suggests another interpretation: it is 
a pasting map between different patches of a larger presymplectic manifold. Let us describe 
an example of such a construction and see, whether it can be of any use or not. 

We shall denote the two maximal extensions of $ defined above with Aq = or Ag~ = 
by $_ and $ + , respectively; their domains in Ti are T>- and T> + f\. Consider two copies of 



Observe that the maximal extensions can be chosen differently, overlapping again at the bound 
trajectories, but the new domain T>'_ containing points of all scattering trajectories expanding to 
the right or contracting from the left and the new T>' + containing points of all scattering trajectories 
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(T E , e E ) denoted by (T' E , G' E ) and (T" E , &' E ), and one copy of (r\, 6i); let $' + : r x ^ T' E has 
the domain V + C r x and let : Ti i— > T'^ has the domain X>_ C lY Let us paste Ti and 

together along V + and Q' + (V + ) by $' + and similarly r x and along V_ and 
by $" . The result is a well-defined (possibly non-Hausdorff) presymplectic manifold, which 
will be denoted by (T,0). 

Problems arise, however, if one tries to find a physical interpretation of this construction 
and considers the role and position of observers. An idea, which seems reasonable, is that 
each dynamical system like Ye or I\ is to describe, in an idealized manner, what a family 
of communicating observers can do with the shell. They can send it contracting in different 
dynamical states, observe its motion and also observe different states of expanding shells. 
It seems further reasonable that such observers could be placed somewhere in the right 
asymptotic regions of A4 + and M._ (or left, but not both right and left, because right 
observers cannot communicate with the left ones). Indeed, the observers are in _ before 
they throw in a shell contracting from right, and then they are in M.+; analogously for shells 
expanding to right, the observers are first in M. + and, after the shell passes, in M.-. In this 
sense, the right asymptotic regions of all spacetime shell solutions of one system are to be 
considered as identical. 

Next consider the pasting of three systems T' E , Ti and Y" E . Nothing seems to hinder us 
in making the assumption that each of the three dynamical systems has its own observer 
family of the above kind before the pasting. It is, therefore, conceivable that some of the 
three families remain distinct after the pasting. Then the observers of T' E might send in a 
shell and this shell disappears behind a horizon for these observers, but appears during its 
motion somewhere else, where it can be observed by another of the three families. Let us 
look to see what happens with the families if we perform the pasting. 

It is a pasting of phase spaces; but such a pasting implies also that points of spacetimes 



that expand to the left or that contract from the right. 
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will be identified. Indeed, let us chose any point (u'_,v'_,P+) G T' E that lies at a bound 
trajectory. (u'_ , v'_ , PL) determines via Eqs. (p6[) and (|39|) the Schwarzschild mass of the 
right subspacetime M' + , via Eqs. fl3~7|), (p8|) , (g3|) and (^5j) a point (£/|,K[_) of where 
the shell is, as well as the four-velocity of the shell at (U' + , K[) represented by P'^ through 
Eq. (PH|) . By 3>!j7 , the point is identified with (u,v,P) G I\ and this, in turn 

by $" with the point (u'L,v'L, P") G T^. Again, (u'!_,v'i,P") determines the mass E'[ = E' + 
of M% the point (U'_(_, V+) of .M'j and a four-velocity at (U'_l, V%). 

It seems, therefore, that the points (U' + ,V+) and (U",V") must be identical after the 
pasting, because they lie at the trajectory of one and the same shell. In fact, one can 
easily see that all points of the trajectory then lie in these particular spacetimes Ai' + and 
Ai'l_ and must be identical in pairs. Thus, and have a one-dimensional set of 
points in common. This is very strange unless the spacetimes and A4'L are themselves 
identical. In fact, there is only one isometry mapping onto M"^ and preserving the 
trajectory points. It is natural to assume that such an identification is performed in all 
cases where the above construction is viable; that is, M' + and + are identical for all 
E' + = E+ G (E- — M,E- + M). Notice that there is no analogous argument for M!_ and 
A4'L, because a point of Y' E (or Y" E ) does not determine a unique point of M.' + (or -M" ). 

It seems to follow that the pasting leads to identification of all three families of observers. 
Indeed, if the observers of the system T' E throw in a shell in a bound state (remember that 
there are bound states passing through any radius, and, moreover, 'lower' bound shells can 
in any case be arraged indirectly by the asymptotic observers) the same must be done at 
the same time and radius by each of the other families. This is 'disappointing' but not 
desastrous. However, if the observers of the system T' E throw in a shell in a scattering state, 
then the same must be done by the observers of Ti, but it must not be done by the observers 
of T E , because contracting scattering states of Y' E are distinct from those of Y" E according 
to the pasting procedure. This seem to be a paradox; I have not found any way out of it as 
yet. Thus, the pasting does not seem to work. 



39 



To summarize: The map $ that realizes the equivalence between the systems (Ti,0i) 
and (Te, ®e) is afflicted with two problems. On one hand, a map satisfying all requirements 
1-3 at the beginning of Sec. [V] does not exist, at least for the case of dust. One might be 
able to weaken the continuity in some cautious way and still preserve the physical content, 
but this is only a speculation that must be studied. On the other hand, the map $ is not 
unique; this does not seem to lead to any serious difficulty as long as just the classical theory 
is concerned, but it can be a handicap for the self-adjoint extension methods in the quantum 
theory. 
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APPENDIX 

In this Appendix, we shall sketch the main steps of the line of reasoning in terms of 
Schwarzschild coordinates so that the important special case E_ = of flat shell interior 
can be dealth with. 

The Warsaw approach starts in this case from the action of the form 

Si = j dt ( PT f + prR - MC X ), (69) 
where T and R are coordinates on a two-dimensional Minkowski half-spacetime R > and 

In order to reduce this action to the Cartan form, we introduce the Schwarzschild (here: 
Minkowski) momentum P$ by 

Ps := -Rarctanh-^. (71) 

The constraint C\ = is then identically satisfied by 
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p T = -M(R)oosh^ + M2{R) 



R 2R 



p R = M(R) sinh^, 



and the Cartan form becomes 



0! = -M(R) ( cosh§ - j rfT + M(P) sinh^-dR. (72) 

\ i? 2ix / i? 



The Potsdam approach in the Schwarzschild coordinates (Ref. [H] ) is awkward, because 
there are 16 disjoint ranges of validity of these coordinates (4 quadrants for each Kruskal 
subspacetime) . We shall need four sign functions a e and b e , e = ±1 and an abbreviation 
sh a x in order that we can catch all 16 quadrants by a single formula. The definitions are: 
a e := sgnF e , where F t := 1 — 2E e /R, b t := +1 in the past and b e := —1 in the future of the 
Kruskal event horizon in A4 e and 

. e x + ae~ x 
sh a x := . 

a 2 

The momentum PJ conjugate to R t is determined by the 'Schwarzschild velocity' dR/dT of 
the shell as follows (cf. Ill): 



P|:=i? e arctanh(l^|) \ (73) 

where dR e /dT e is the derivative of the Schwarzschild radius R e with respect to the 
Schwarzschild time T e along the shell from the e-side. The action reads 

S 2 = Jdt {[P S R - Ef] + V{R + - R_) - vC 2 ) , (74) 

where v and v are Lagrange multipliers and 



Ps 

bRJ F sh, 



C "' : R 

There is a secondary constraint, x = 0) where 



+ M(R). (75) 



x ■-- 
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explicitly, 



X 



bJ\F\ sh_ a 



Pa 
R 



The three constraints C2 = 0, R+ = R- and x = can De rewritten as the following 
system of equations for R + , E + , P + a + and b + : 

R + = R, a + = sgnF + , 



■a+ 



- smh 



R 



R ' 



we have left out the index "— " and we have already set E = 0, a_ 
F- = 1 everywhere. It follows immediately that 



cosh- 



M{R)\' 



R R 



sinh 2 ^, 
R 



-1, 6_ = +1 and 



(76) 



and 



tanh-^- 
R 



Eq. flT6"| ) can be rewritten in the form 



sinh-^- 



a+ 



M(R) 



. cosh-|- R 



(77) 



F+ — Q1Q2, 



(78) 



where 



1 _ M{Rs) p-/ R 
R ' 



Q 2 



(79) 
(80) 



are two important abbreviations. We also define the following sign functions: 



qi : = sgnQi, g 2 := sgnQ 2 - 
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Then, Eq. ( |T8|) implies 



a+ = <7i<?2, 



(81) 



and 



J5+ = |(i-QiQ a ). 



Explicitly, 



E + = E_ + M(fl)coshf _ JOS 



(82) 



Let us turn to Eq. (f77|) . We use the following simple identity 



arctanhr a = - In fa — — — 
2 V 1-x 



that holds for a = ±1 and all \x\ < 1. It implies that 



However, 



a + 



Qi 



Qi 



Q2 



because of Eq. (^T|) so that 



P+ = -In 



Qi 



or equivalently 



\ 



Qi 



e Ps/R 



(83) 



As for b + , we just substitute Eq. (|83|) for e ±p s / R and Eq. (f7§D for F + into the equation 
X = 0; this leads to 



6+ (|Qi|e p * /R - a+|Q 2 |e- p s /«) = 2sinh^ 



i? 
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Using the definitions of qi, q2, Qi and Q2, we obtain easily 



(84) 



From Eqs. (|81|) and fl34|), it follows that Qi changes sign at the "black hole" and Q2 at the 
"white hole" horizon of 

Subsituting Eq. (|82"D for E + , Eq. (|83D for Pg and R for R + and i?_ into the variational 
principle we obtain the reduced action with the Cartan form 0^ given by 

Qi 



2 := - In 

2 



Q2 



dR - M(R) ( cosh% - ^P-] dT + . 



(85) 



R 2R 

The cyclic coordinate T_ is automatically excluded, because it is contained only in one term 
of the form E_dT_ and E_ = 0. The manifold is determined by the inequalities 



R > 0, R ^ 2£ + , E + > 0, 

that must be satisfied by the coordinates R, Pg, and T + ; it consists of four open disjoint 
submanifolds. 

Finally, we have to show that there is a transformation 

T + = 4>(T,R,P S ), 



such that Qe becomes to 0i plus possibly a closed form, if we substitute <fi for T + . In an 
analogous way as in Sec. |V A|, we obtain 



4>(T,R,P S ) = T + A(R,P S 



and 

dA 2E+R 2 - M 2 (R)R 2E+R + M 2 (R) 
— = sgnP 5 === - sgnP s == . 

dR 2(R-2E + )^V{R) yV{R) 

In the case of dust, M(R) = const, the differential equation (|86|) reduces to 

dA 2El - m 2 /•;. i 1/; 2 - M 2 ) 

— = sgnP s ± + sgnP s ■ 

dR JV(R) 



(R-2E + )JP(R) 



(86) 



(87) 
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where 

V{R) := (El - M 2 )R 2 + M 2 E + R + M 4 /4. 

The integration of Eq. (|8T| ) is completely analogous to that of (|61| ) and the results are 
also analogous, only the divergence of A 2 at R = is shifted to R = 2E + . This divergence 
of A2 at R = 2E + does not lead to any problem; it originates from the singularity of the 
coordinate T + on the Kruskal spacetime Ai+ rather that from any geometrical effect. 
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